Recombination fluorescence in ultracold neutral plasmas 
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We present the first measurements and simulations of recombination fluorescence in ultracold 
neutral plasmas. In contrast with previous work, experiment and simulation are in significant dis- 
agreement. Comparison with a recombination model suggests that the disagreement could be due 
to the high energy portion of the electron energy distribution or to large energy changes in elec- 
tron/Rydberg scattering. Recombination fluorescence opens a new diagnostic window in ultracold 
plasmas because it probes the deeply-bound Rydberg levels, which depend critically on electron 
energetics. 
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Ultracold neutral plasmas can provide new insights 
into relaxation phenomena in strongly coupled Coulomb 
systems. They are created by photoionizing laser cooled 
atoms and therefore have very precisely controlled initial 
temperatures and densities [H,13, Photoionization 
causes an impulsive hardening of the interparticle po- 
tential and the pathway to (non)equilibrium can be 
studied in detail @, @] • 

The electron system equilibrates on the shortest time 
scales and lar gely determines how the plasma expands 
[1> S E3) II > K2] • Indirect measurements of the elec- 
tron temperature agree well with simulations 0, Hjj , 
Ibi ] . Disorder-induced heating, three body recombination 
(TBR), and electron/Rydberg scattering are important 
heating mechanisms in the plasma. The latter two are 
thought to occur on times that are long compared to 
the electron plasma frequency. Thus TBR can probe the 
electron system at early times in the plasma evolution 
after disorder- induced heating has occurred [l5j |. 

In higher temperature plasmas, the TBR picture is an 
electron and ion colliding in the presence of another elec- 
tron. The TRB rate is the two-body collision frequency 
(navth) multiplied by the probability that a third body is 
one collision distance away (nb 3 ). The collision distance 
and cross section depend on the electron temperature 
(b = e 2 1 iTTeoksT and a = Trb 2 ), giving the well-known 
TBR rate 03 = Kn 2 T~ 9 / 2 , where TZ is the TBR rate 
coefficient [161 ]. 

In a strongly-coupled system, this binary collision 
picture breaks down. When the nearest-neighbor 
potential energy exceeds the kinetic energy (r = 
e 2 n 1 / 3 /47reofc_BT > 1), the plasma is strongly coupled. 
The average distance between electrons becomes compa- 
rable to the collision distance b ~ n" 1 / 3 . The electrons 
are in a constant collision and TBR becomes, in a sense, 
many-body recombination. A numerical simulation of 
early-time recombination suggest that even moderate ini- 
tial coupling in the electron system changes the recom- 



bination rate by perhaps a factor of two [17fl. Other the- 
oretical treatments suggest larger changesjlil. Hij . How- 
ever, there is no experimental evidence that TBR departs 
from the standard formulas (ioj . 

In this paper we present a new study of three-body 
recombination in ultracold neutral plasmas. We measure 
the time-dependent fluorescence emitted by the plasma 
following electron/ion recombination for a range of ini- 
tial plasma densities and electron temperatures. The flu- 
orescence signal is the end result TBR, electron-Rydberg 
collision, and radiative cascade. Thus, it is sensitive to 
the details of the evolving electron system. We compare 
these results with a numerical simulation and find signifi- 
cant disagreement. A recombination model suggests that 
this is related to the high-energy portion of the electron 
energy distribution or to large energy changes in elec- 
tron/Rydberg scattering. 

Our ultracold plasmas are created by photoionizing 
laser-cooled calcium atoms in a magneto-optical trap 
[H HH. The plasma density profile is approximately 
Gaussian and spherically symmetric with the density 
n(r,t) = no exp(— r 2 /w 2 (t)) and w(0) — 0.5 mm. The 
initial ion temperature is equal to the neutral atom tem- 
perature (~1 mK). The initial electron energy is equal to 
the difference between the photoionization laser photon 
energy and the atomic ionization energy. The lower limit 
is set by the laser bandwidth and by continuum lowering 
effects to approximately 2E e /3kB = 0.5K. Our two-step 
photo-ionization process ionizes nearly 100% of the atoms 
in the MOT, and our maximum plasma density is set by 
the MOT density to n < 2 x I0 10 cm" 3 . 

Fluorescence from the plasma is collected by a fast lens 
system and detected using a PMT. An interference filter 
blocks light from the lasers and the MOT. Most of the 
fluorescence results from emission at 410 nm. A typical 
fluorescence signal and partial energy level diagram are 
shown in Fig. Q] At early times, the fluorescence sig- 
nal grows as recombination proceeds. After a few fj,s the 
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FIG. 1: Top panel: Partial energy level diagram for neutral 
calcium. For a fraction of recombination events, the radia- 
tive cascade produces fluorescence at 410 nm. Lower panels: 
Recombination fluorescence signal at 410 nm for T = 57K 
and three different densities, experiment (black points, left 
scale) and simulation (gray line, right scale). The simulation 
changes more much rapidly with density than the experiment. 



plasma expands and the density decreases. The fluores- 
cence signal falls as recombination slows down and the 
plasma moves out of the view of the fluorescence collec- 
tion optics. 

The fluorescence signal is the end product of a se- 
ries of collision events 



22l | . Free plasma electrons and 



recombined Rydberg atoms are in a quasi-equilibrium. 
Collisions form and then ionize Rydberg atoms, most of 
which have a binding energy roughly equal to fc^T. This 
collisional equilibrium gives rise the so-called "thermal 
bottleneck," which is a minimum in the Rydberg state 
distribution near — 4fc^T 17(. A collision occasionally 



occurs that puts the atom in a Rydberg state below the 
bottleneck. Additional collisions and radiative cascades 
move a fraction of these atoms into the fluorescing state. 

We have performed a numerical simulation of recom- 
bination and fluorescence in the plasma using a method 
similar to Refs. |23j|. Like the experiment, we assume 
the plasma to be nearly charge neutral with a spheri- 
cally symmetric Gaussian density distribution. The elec- 
trons are in thermal equilibrium at a temperature T. 



The plasma expands self-similarly, and the density can 
be written as n(r,t) = no exp[— r 2 /w 2 (t)]. The equations 
of motion at this level of approximation are described 
in Sec. IIIc of Ref. 23]. The main trends are that the 
expansion decreases the thermal energy of the electrons 
while increasing the radial speed of the ions. When re- 
combination occurs, the density decreases slightly and 
the temperature increases. An electron collision with an 
atom that decreases the atom's energy increases the elec- 
tron temperature and vice versa. 

We compute the evolution of the atom population us- 
ing a Monte Carlo technique. During a time step, we 
compute the probability of formation of an atom at each 
position r and which state it is in using the method 
of Sec. IV of Ref. 23j. The method gives more atoms 



formed at the center of the plasma where the density 
is high and the atoms are mostly formed in high Ryd- 
berg states 2J, |25j. The radial velocity of the atom is 
taken to be the radial velocity of an ion at the position 
the atom is formed; after the atom is formed, we track 
its position and velocity so we know the correct electron 
density near each atom and we know whether the atom 
is within the detection region when the photon is emit- 
ted. After the atom is formed, its state changes due to 
electron collisions and radiative decay. The electron col- 
lisions can reionize the atom, excite it to higher energies, 
or de-excite it. The electron collision processes are com- 
puted with a Monte Carlo technique using the rates of 
Ref. I2al. One difference with Ref. [231 is that we used 



the exact radiative decay rates for nl 



'£' for n < 30. 



We did this because the £ states are not necessarily sta- 
tistically mixed for lower n and the fluorescence signal 
depends on getting the cascade correct. When an atom's 
principle quantum number drops below 5 it is counted as 
a fluorescence photon and removed from the simulation. 
The calculation has only one free parameter related to 
the location of the thermal bottleneck. The formulas in 
Refs. (Him lead to a bottleneck at -2k B T. Changing 
this value changes the shape of the fluorescence curves 
somewhat because the details of the fluorescence signal 
depend on the high energy portion of the electron distri- 
bution. However, there was no choice for the bottleneck 
that allowed us to reproduce all of the measurements. 

The simulation and experiment are compared in the 
lower three panels of Fig. [1] In these data, the initial 
electron temperature is T = 2i? e /3fcs = 57 K, and the 
initial plasma size is w(0) = 0.5 mm. Qualitatively, the 
simulation and experiment are in reasonable agreement. 
For the lower two densities, the rising edge of the simu- 
lation and experiment have the same time dependence. 
However, there are quantitative mismatches. The rela- 
tive change in signal strength with density is completely 
different for the simulation and experiment. In addition, 
the time dependence of the high density simulation does 
not match the experiment. 

This disagreement is surprising. Previous work has 
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demonstrated good agreement between this simulation 
and the experiment. This simulation and related ap- 
proaches give the proper plasma expansion velocity, elec- 
tron temperature, and even the high n Rydberg state dis- 
tribution. However, recombination fluorescence opens a 
new window in these ultracold plasmas because it probes 
the deeply-bound Rydberg levels. The population in 
these levels depends critically on the high-energy portion 
of the electron distribution. 

Some insight can be gained by comparing to a simple 
rate equation model. The rate at which Rydberg atoms 
are formed in the plasma can be written as 



~db 



= -S 



(1) 



where S is the collisioinal ionization rate, 1Z is the three- 
body recombination rate coefficient, n r (n e ) is the Ryd- 
berg atom (electron) density, and T is the electron tem- 
perature. This simplified model treats all Rydberg levels 
as essentially equal and neglects transitions between dif- 
ferent levels. It neglects radiative decay and assumes the 
plasma is charge neutral. Recombination increases the 
electron temperature, and the corresponding rate equa- 
tion can be written as 



die 

dt 



1 dn r 
n e dt 



T e + TqTq 



n e (0) 



1/3" 



(2) 



where Tq, Tq, and n e (0) are the initial electron temper- 
ature, strong coupling parameter, and electron density. 
The total energy gained in one recombination event is 
the sum of the colliding electron energy and the correla- 
tion energy. We use the relation n r + n e — n e (0) to write 
these coupled differential equations for n r (t) and T e (t) in 
terms of known parameters (27j . 

The relationship between the Rydberg atom density 
and the fluorescence emission rate is non-trivial. How- 
ever, there is some evidence from our simulation that a 
single collision de-excites the Rydberg atom to radiating 
states. This suggests that the Rydberg atom density at 
a time t\ is proportional to the fluorescence signal s(t) 
integrated to a time £2, 



n r (ti 



s(t)dt. 



(3) 



A comparison of this simplified rate equation model to 
the experiment is shown in Fig. [5] For this comparison, 
we choose ti early enough that the plasma has not started 
to expand. 

The relatively good agreement between the rate equa- 
tion model and the data and the poorer agreement of 
the simulation deserve some attention. The simulation is 
much more complete than the model. It uses the most 




FIG. 2: A plot of Eq. [3] comparing the rate equation model 
(lines) to the experiment (points) for different temperatures 
and densities using ti = ti = 1.5 /is and a single overall 
proportionality constant. 



widely accepted approach for Rydberg atom formation, 
for Rydberg/electron collisions, and for Rydberg state 
distribution. It uses the exact hydrogenic decay rates. It 
also includes geometrical effects of the fluorescence collec- 
tion optics. On the other hand, the rate equation model 
doesn't include any effects due to Rydberg state distri- 
bution or to the cascade down into fluorescing states. 
Perhaps, then, the simulation tells us something about 
effects beyond the "standard" cold plasma model. It is 
likely that these are related to details of the high-energy 
portion of the electron distribution. A non-Maxwellian 
tail would change the location of the thermal bottleneck, 
the deeply-bound Rydberg state distribution, and the 
cascade rates. 

A simulation in Ref. [17| suggested that strong cou- 
pling in the initial electron system might change the 
three-body recombination rate by perhaps a factor of 
two. Guided by the rate equation model, we can ex- 
plore this possibility. For early enough times, the col- 
lisional ionization term, S, in Eq. [T] can be neglected 
because the recombined Rydberg atom density is small. 
With this assumption, the Rydberg atom density is 
n r (t)/n e (Q) ~ n 2 e {0)T~ 9 / 2 t. Using Eq. [3] and the defi- 
nition of the strong coupling parameter, we can write 



r 



-3/2 



s(t)dt 



2/9 



(4) 



For the low density, high temperature data, the strong 
coupling parameter can be reliably calculated from 
the initial experimental values 14]. However, when 



the plasma is generated closer to the strongly-coupled 
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bution is different from that predicted by the standard 
cold plasma simulation. We note that when the ultracold 
plasmas start with the electron strong coupling parame- 
ter in the range of 0.1 < T < 10, the early time evolution 
of the electron and Rydberg state distributions are not 
well known. Measuring the plasma fluorescence provides 
an opportunity to study this regime. 

This work is supported in part by Brigham Young Uni- 
versity, the Research Corporation, the National Science 
Foundation (Grant No. PHY-0601699) and the Chem- 
ical Sciences, Geosciences, and Biosciences Division of 
the Office of Basic Energy Sciences, U.S. Department of 
Energy. 



FIG. 3: Strong coupling parameter T from Eq. [4] plotted vs. 
(a) initial calculated F and (b) the temperature-adjusted F. 
The experimental data is taken from Fig. [2] In panel (b) 
the adjusted temperatures are adjusted as follows: Black cir- 
cles (T = 57K) no adjustment, upright triangles (T = 30K) 
adjusted to 38K, inverted triangles (T — 0.5K) adjusted 
to 20K. The solid line plots the equality F (experiment) = 
F (adjusted). 



regime, the electron temperature rapidly increases and 
the density decreases to values that can be significantly 
different from the initial ones. However, it is expected 
that the recombination rates are smooth functions of the 
strong coupling parameter. 

In Fig. GHa) we plot the strong coupling parameter 
extracted from the experiment vs. the calculated value. 
The experimental data is scaled to agree with the calcu- 
lated value for the smallest T. The higher T data are in 
obvious disagreement. However, by adjusting the elec- 
tron temperature in the calculation, the experimental T 
can be made to vary smoothly as shown in Fig[3jb). Be- 
ginning at r w 0.1 the experimental T falls well below 
the adjusted value. This is what one would expect if the 
three-body recombination rate is depressed due to many 
body interactions. Of course, this is also what one would 
expect if the electron temperature is simply higher than 
expected. An independent and accurate electron temper- 
ature measurement at early times is obviously needed, 
although if the tail of the electron distribution is non- 
Maxwellian, the meaning of the temperature measure- 
ment becomes less clear. 

In conclusion, we have demonstrated a new fluores- 
cence technique to study three-body recombination in 
ultracold neutral plasmas. This method probes the evo- 
lution of deeply-bound Rydberg states. In contrast with 
previous work, where the evolution of plasma parameters 
can be predicted using simple fluid models, we find sub- 
stantial disagreements. This work suggests that the high 
energy electron distribution or the Rydberg state distri- 
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